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We estimate the expected behavior of a quantum model of heat engine when we have incomplete 

information about external macroscopic parameters, like magnetic field controlling the intrinsic 

energy scales of the working medium. We explicitly derive the prior probability distribution for 

. . . these unknown parameters, a^, (i = 1, 2). Based on a few simple assumptions, the prior is found to 

04 , be of the form n(ai) oc l/oi. By calculating the expected values of various physical quantities related 

to this engine, we find that the expected behavior of the quantum model exhibits thermodynamic- 

^^ ■ like features. This leads us to a surprising proposal that incomplete information quantified as 

Cn I appropriate prior distribution can lead us to expect classical thermodynamic behavior in quantum 

$_^ , models. 
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I. INTRODUCTION 

^- n .....nn 

O I In a recent paper [1], a connection was observed between the notion of prior probabilities [2|, |3| and the optimal 

"^ ] performance of certain models of quantum heat engines [4] . As a central feature of Bayesian inference, these probabil- 

2 . ities quantify the uncertainty or equivalently, the lack of complete knowledge about some parameters of the problem. 

O^ It was found [ij that a power-law choice of the prior, yields the expected optimal performance of the quantum engine 
at certain well-known thermal efficiencies such as Curzon-Ahlborn (CA) efficiency [5]. This observation is remarkable 

^ ' from at least two viewpoints: i) the CA efficiency is well-known to be a finite-time thermodynamic result, where it 

^: . . . . 

f^ , arises within endoreversible thermodynamic models as efficiency at maximum power. But the cycle considered in Rcf. 

^y-, ■ [l| runs infinitely slow; ii) the derivation of this result based on Bayesian approach suggests a very different origin of 

[■~~. ■ these efficiencies by connecting them with the notion of incomplete information. 

I ■ Let us clarify about the source of the assumed incomplete information. Usually, incomplete information about the 

. . , state of a quantum system is meant to refer to an uncertain preparation procedure, the state being defined by a 

>• . 

• i-H . mixed state density matrix p — ^iPi\ipi){ipi\. This ensemble is characterised by a hamiltonian H{a) of the system, 

^ ■ • 

5h ' ^ function of externally controllable parameters (collectively denoted by a) , such as magnetic field, coupling strength 



with environment/subsystems and so on. One may say that the state or the ensemble is characterised by the choice 
of these external parameters. In the present discussion, we assume an ignorance about the exact values of these 
parameters. Our focus is on the following issue: can someone who is ignorant of the exact values of certain control 
parameters for the system, make some reasonable estimates about its behavior? The particular physical problem that 
we choose is the performance of heat engines. In particular, we ask how should one quantify one's ignorance about 
the configuration of a quantum heat engine and what are the estimates for its expected performance? 

Now we may possess some prior information about a parameter which is uncertain in the above sense, such as the 
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possible range of values or, other inputs from the physics of the problem. Based on that we assign probabilities for 
the likely values of this parameter. However, these probabilities are to be interpreted in the sense of rational degree 
of belief [3]- Thus essentially we interpret our ignorance about parameters of the system from a Bayesian perspective. 
Formally, one may say that we are considering an ensemble of subensembles, where each subensemble is characterised 
by a hamiltonian H{a), and the possible values of a over different subsensembles are described by a prior distribution. 
We may refer to the ensemble so considered as the prior ensemble, but the ensemble is regarded here as a theoretical 
construct only. Our aim is to quantify the prior information about the uncertain parameters as a prior distribution 
and to use this prior to make estimates about the behavior of the system. 

The purpose of this paper is to derive and justify the choice of the prior from more basic underlying assumptions 
and to further clarify the subjective approach. Our quantum model of a heat engine consists of a pair of weakly 
interacting two-level systems, which is closely analogous to the model of [l|, |j]. There are two reservoirs characterised 
by temperatures Ti and T2. Additionally, we have two energy scales, denoted by ai{i = 1,2) refering to the level 
spacing of a two-level system. These spacings are controllable externally e.g. via magnetic field, but are here assumed 
to be uncertain. It is these energy scales for which we assign a prior. The form of the prior we derive for parameter 
ai, is the l/a^ law. We show the consequences of using such a prior for the expected values of physical quantities 
such as efficiency, work and so on. We observe that the assignment of the prior is such that it leads us to expect a 
thermodynamic like behavior in these quantum models. 

The paper is organised as follows. In Section II, we present the model for heat engine and summarise its main 
features. In Section III, we assume incomplete information about internal energy scales of the working medium and so 
derive the prior for them based on some simple assumptions. In Section IV, we apply the prior so derived to estimate 
thermal efficiency of the engine, final temperature, work extracted and so on. In Section V, we highlight a specific 
asymptotic limit in which the expected behavior becomes especially simple. Section VI, addresses a special case by 
including a constraint. The final Section VII is devoted to conclusions and future outlook. 

II. QUANTUM MODEL FOR WORK EXTRACTION 

Consider a pair of two- level systems labeled R and S, with hamiltonians Hn and Hs having energy eigenvalues 
(0, ai) and (0, 02), respectively. The hamiltonian of the composite system is given hy H = Hn iS I + I ® Hs- The 
initial state is pini ~ pR® Ps, where pE. and ps are thermal states corresponding to temperatures Ti and T2 (< Ti), 
respectively. Let (ri, r2) and (si, S2) be the occupation probabilities of each system, where 

with ^2 = (1 — ri) and S2 = (1 ^ si)- We have set Boltzmann's constant fee = 1. The initial mean energy of each 
system is 

where z = 1,2 denote system R and S respectively. Within the approach based on quantum thermodynamics [6|- 
8|, the process of maximum work extraction is identified as a quantum unitary process on the thermally isolated 



composite system. It preserves not just the magnitude of von Neumann entropy of the composite system, but also 
all eigenvalues of its density matrix. It has been shown in these works that for ai > 02, such a process minimises the 
final energy if the final state is given by pun = Ps ® PR- Effectively, it means that in the final state the two systems 
swap between themselves their initial probability distributions. The final energy of each system at the end of work 
extracting transformation is 

where i ^ j. The average work per cycle defined as W = Tr[(pini — PRn)H] = Eini — Efin, is given by 

1 1 



W{ai,a2) = (ai - a2) 



(l + e'^i/^i) (l + e°=/^2) 



(4) 



To complete the cycle, the two systems are brought again in thermal contact with their respective reservoirs. In this 
cycle, the heat extracted from the hot reservoir is 



h = fli 



1 1 



(l + e^^/Ti) (i + e°2/T2) 



(5) 



The efficiency of this engine 77 — W/Qi is 



77=1-—. (6) 



Note that for 02 = ai,W — 0,Qi > and 77 = 0; for 02 ~ ai{T2/Ti), we have the limiting values of W^ = and 
Qi = and 77 = I — {T2/T1). The operation of the machine as a heat engine defined asW>0 and Qi > 0, is satisfied 
if 

ai{T2/Ti) <a2<ai. (7) 

III. BAYESIAN APPROACH 

Now consider a situation in which the temperatures of the reservoirs are given a priori such that Ti > T2, but about 
the parameters Oi and 02 we only know that, 

• Oi and a2 represent the same physical quantity, which is level spacing for system R and S respectively, and so 
they can only take positive real values. 

• If the set-up of R+S has to work as an engine, then criterion in Eq. ([T]) must hold, whereby if one parameter is 
specified, then it redefines the range of the other parameter. 

Apart from the above conditions, we assume to have no information about ai and 02- The question we address in 
the following is: What can we then infer about the expected behaviour of physical quantities for this heat engine 
(such as work per cycle, efficiency and so on) ? We shall follow a subjective approach to probability to address this 
question. This implies that an uncertain parameter is assigned a prior distribution, which quantifies our priliminary 
expectation about the parameter to take a certain value. We denote the prior distribution function for our problem 



by 11(01,02). The prior should be assigned by taking into account any prior knowledge or information we possess 
about the parameters. For example, if ai is specified, then the prior distribution for 02, 7r(a2|ai) is conditioned on 
the specified value of Oi, and is defined in the range [ai9,ai\, where = T2/T1, because we know the set-up works 
like an engine if we implement Eq. ([7]). 

The expected value of any physical quantity X which may be function of oi and 02, is defined as follows: 

X= XIi{ai,a2)daida2. (8) 



A. Assignment of the prior 

In Bayesian probability theory, the assignment of a unique prior is a central issue. It should quantify not only the 
prior knowledge about the parameter in the particular context, but should meet a consistency criterion according to 
which different observers in possession of equivalent information should assign similar priors. Suppose, one knows 
only the range within which the parameter takes its values, then intuition suggests that a uniform distribution may 
reflect the state of our knowledge. But such a choice is not invariant under reparameterizations. Our aim in the 
following is to motivate the assignment of prior distributions for level spacings ai and 02, when the physical problem 
at hand is a heat engine as described in Section plj) . 

It seems convenient to speak in terms of the two observers A and B, who wish to assign priors for ai and 02- The 
assignment is based on the following assumptions: 

(a) Same functional form of the prior is assigned to ai and a2 in the initial state, denoted by n(ai) and 11(02) 

respectively. Further, we assume that the prior can be expressed as Tl{ai) oc df{ai)/dai, using a continuous 
differentiable function /(a,) {i — 1,2). 

(b) The conditional prior distributions p{aj\ai), implying distribution for aj given a value of a.i, or p{ai\aj) for the 

converse case, have the same functional form as above, and we assume that p{aj\ai) oc df{aj)/daj, where the 
dependence on ai may be present in the normalisation factor. Similarly, we assume p{ai\aj) oc df{ai)/dai. 

Assumption (a) is reasonable since both ai and 02 represent the same physical quantity, and the state of knowledge 
of A and B about them is the same, which is the fact that their values in the initial state lie in a preassigned range 
[oniin, Omax]. For simplicity and symmetry, we take this range to be identical for ai and 02. Presumably, this range 
depends on the experimental setup, and we assume similar apparatus for controlling the level spacings of systems R 
and S. 

Now if one of the parameters is specified to an observer, say ai to A, then A knows that the machine works as an 
engine only if the range of 02 is [aiO, ai], where ai here represents some fixed value. Even so, assumption (b) states 
that A must assign the same functional form to the prior for 02 as it used for assigning to ai. Although this does 
not represent the general case, we assume this for simplicity and in the following analyse the consequences of these 
assumptions. 

Thus from (a), we have 

1 df{a^) 



where the normalization constant is determined as 

M = 



df{ai 



■dai 



dcLi 



Following assumption (b), the conditional probability distributions are given by 

1 df{a2) 



and 



where 



n(a2|ai) 



n(ai|a2) 



Ni 



Ni da2 



1 dfjai 
N2 dai 



(10) 



(11) 



(12) 



dfia2] 



da2 



= /(oi) ~ f{ait 



(13) 



and 



N-2 = 



02/9 



d/(ai) 



dai 



dai 

f{a2/e)-f{a2), 



(14) 



are the respective normalization constants. Now using the product law of probabilities, the joint prior 11(01,02) as 
expressed by observer A 



11(01,02) 



n(o2|oi) • n(oi) 

df{a2)/da2 

[/(ai) - f{aiO)] [/(Omax) - /(flmin)] ' 



df{ai)/dai 



(15) 



or equivalently in terms of B, 



11(01,02) 



n(ai|o2) • 11(02) 

df{ai)/dai 



d/(02)/d02 



(16) 



[f{a2/0) - /(02)] [/(Omax) - /(flmin)] ' 

As shown in the Fig. 1, each observer assigns different ranges of values to oi and 02. Each can use its own joint prior 
to make estimates about any quantity. Now since either approach, that of A or B, is equivalent, consistency would 
require that they make similar estimates as to a given quantity. Also it is reasonable to assume that for some given 
pair of values (01,02), which are in the allowed range of each observer, each of them should assign same probability 
for choosing oi within the small interval doi around the given value oi, as well as of choosing 02 within the small 
interval do2, around the given value 02. Clearly, for such a pair, both oi and 02 have to lie in the interval [omin, 
Omax]- One such pair of values, which is definitely common to both A and B would be (ai,oi) i.e. when 02 = oi. 
Thus equating the probabilities assigned by A and B for this case, we obtain from Eqs. ()15|) and ([T6l) . 



/(oi)-/(oi0) = /(oi/e)-/(oi). 



(17) 
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FIG. 1: Observer A (B) assigns the same range of values for ai (02) in the initial state of spin R (S). But the range assigned 
for the other parameter a2 (ai), conditional on the operation as an engine, is different. This fact manifests such that A and B 
in general arrive at different estimates for physical quantities. 



which can be rewritten as 



2/(ai)-/(ai0) + /(ai/0). 



(18) 



This functional equation has a unique solution, f{x) — Alnx, upto an additive constant. Thus we find the explicit 
form of prior (Eq. ([5])) as 

1 1 



n(a.) 



(19) 



V am,„ ) 

which is functionally the same as Jeffreys' choice j2[, also employed in a previous study [l|. The joint prior is then 
given by 

1 1 



11(01,02) 



In 



(i)i^fe)^"^"^)" 



(20) 



The joint prior derived above is relevant only with regard to the final states of R and S. As discussed in the previous 
section, the initial state of system R (S) depends only on parameter ai (02). This fact will be used in the following 
in calculations on the expected values of quantities. 

IV. EXPECTED VALUES OF PHYSICAL QUANTITIES 



In this section, we use the priors assigned above, to find expected values for various physical quantities related to 
the engine. For this purpose, we employ the definition in Eq. ([8]). These expected values refiect the estimates by an 
observer who is assigning the priors. As observed in Eqs. ()15p and ()16|) . there are two methods of writing the joint 
prior. So in principle, there are two ways to calculate the expected value of some quantity and its value will depend, 
in general, on the method used. 



A. Internal energy 



We calculate the expected values of internal energies for systems R and S. These values can then be used to find 
the expected work per cycle, heat exchanged and so on. 

(i) Initial state: For a given a^, the internal energy E^l^- is given by ^. The expected initial energy is defined as 






(21) 



where i = 1,2. Note that _Ej„j depends only on a^, so we need to average over the prior for a.^ only. Thus we obtain 

-1 



or explicitly 



E. 



{i) 



E 



In 



■(i) 



In 



On 



dai 



1 + e'^„i„/T. 



(^^inax Q-ininj \ -Li lH 



(22) 



(23) 



^ 1 + e''""-/^' 

(ii) Final state : In this case, the internal energy of R as well as S, is function of both ai and 02 (see ^) and so 
the expected values are obtained by averaging over the joint prior, n(ai, 02). For instance, the expected final energy 
of system S (denoted by superscript (2)) as calculated by A, 

1 



4™(A) - K 



—— — I — dax / da-i. 



K(\ 



(Omax - Oniin) + Tl In 



1 + e°"""/^^ 
1 + e'^"--/^! 



Similarly if calculated by B, 



E%m = K 



da-2 



ailB 



da I 



(24) 



(25) 



,a, (l + e^'i/^Oai' 
where K = [In (1/0) ln(amax/amin)]~^- The latter integral cannot be completely solved. To simplify, we rewrite it as 



EfL{B)^K 



dao 



1 



(26) 



„, (l + e°i/^i)ai 

Considering the inner integral as the first function and unity as the second function and integrating by parts, leads to 

r^/^ dai 

'02 



E'fi(B) =. Ka, 



-K 



(1 + e'^i/'^i)ai 



02 



dai 



da2Ja2 (l + Ci/'^Oai 



da2- 



Here we use Leibniz integral rule 



— / f{x)dx 

dy Jg(y) 



—, — nKv)) — -7—na{y)), 

dy dy 



(27) 



(28) 



to solve the second term of the Eq. ([27|) and to finally obtain 



dai 



—(2) /"'"/'' 

E)l{B) ^ Ka2 ■ 

•Jao 



-K 



To In 



(1 + e°i/'^i)ai 

1 + e""""/'^2 ■ 
1 + e'^--/^^ 



Tiln 



1 + e°-°/^^ 
1 + ea„,,/ri 



(29) 
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In general, the expected final energies of S, as given by Eqs. (l24l) and (l29l) according to A and B, respectively, are 
not equal. One would expect that if the state of knowledge of A and B is similar, then they should expect the same 
value for a given quantity. However, the difference in the values expected by A and B is not so surprsing in light of 
the fact that different ranges for variables ai and 02 are being employed by them, as shown in Fig. 1. 

Similar feature is also observed in the expressions for expected energy of system R (supersript (1)), which we provide 
below for sake of completeness. 



4m (A) = y y 4™n(a2|ai)n(ai)da2ciai, 



which implies 



4»(A) = K 



da I 



da-i 



It is interesting though to observe that the above integral is identical to the one in Eq. 



5(1 



(2) 



E),UA) = E)/JB) 



The second method however, yields 



^W 



E)L{B) = K 



K 



da2 



ailB 



1 



1 



(1 + £"2/^2)02 7„, 



Tain 



1 



oa„i„/T2 



(30) 



(31) 



(32) 



(33) 



_ 1 + e°".-/T'2 

In the next section, we look at these expressions in a particular limit in which the expected values obtained by the 
two observers yield similar results, so that a meaningful analysis can be carried out in this limit. 



ASYMPTOTIC LIMIT 



As remarked above, observers A and B are supposed to arrive at similar conclusions. So they should arrive at 
similar estimates for physical quantities using their respective priors. This happens in the limit, when Omin << T2 
and Cmax >> Ti. In this limit, Eq. (P^ is approximated as 



E. 



■(«) 



In 2 



(34) 



ln(2m^) '■ 

^ dmin ' 

The ratio (amax/amin) in the above may be large in magnitude, but is assumed to be finite. 

Similarly, it is remarkable to note that in this limit, not only the expected energy of a system (R or S) calculated 
by either of the methods (A or B), is the same but also its value for system R or S is also equal. In particular, the 
first term of Eq. (j29|) can be shown to be negligible in this limit. Thus we have (omitting the observer index) 

,W In 2 (l-0)Ti 



E\ 



(35) 



''™^M!— ) ln(i) ' 

where i — 1,2. Further insight into this may be obtained if we estimate the final temperatures of systems R and S 
after the work extraction process. Now if values of both a\ and a-i are specified, the temperatures (T/) of the two 
systems after work extraction, are given by 



T[ = T, 



02 



and Ti = Ti 



02 
ai 



(36) 



In general, the two final temperatures are different from each other. Within the present framework, when we look at 
the expected values of the final temperatures as calculated by A or B, we find 

(i-e) 



Ti=T^ = Ti- 



(37) 



ln(l/6l)' 

It is interesting to find that the assignment of the prior is such that the two systems are expected to finally arrive at 
a common temperature. Going back to Eqs. (1341) and (1351) for the energies, we see that they satisfy a simple relation 

i?j„j oc Ti and -Bfj„ oc Tj. This is reminiscient of the thermodynamic behavior of a classical ideal gas. 

7=r ^(*) -T^(^) 7T 



Next, the heat exchanged between system i and the corresponding reservoir is given by Qj = E. 



Efm- Qi > 



(Qj < 0) represents heat absorbed (released) by the system. Then the expressions for the heat exchanged with the 
reservoirs in the said limit are as follows: 



In 2 



In 



and 



Q2 






In 2 



1 



In 



ln6l 



6* In 6* 



Ti, 



(38) 



T, 



(39) 



Now the expected work per cycle is defined as: W = Qi + Q2- Thus the efficiency may be defined as 77 = 1 + Q2/Q1 
Explicitly, using Eqs. dSH]) and dSH]) we get 

6'ln6i+(l-6i) 



?7= 1 + 



(40) 



ln6i+(l 

This is the efficiency at which the engine is expected to operate for a given 9. The above value is function only of the 
ratio of the reservoir temperatures and is plotted in Fig. 2. In the limit of small temperature differences, 

(1-61) (l-^^2 



77 



+ 0(1 -61)^ 



(41) 



3 9 

Before closing this section, we note that the constant of proportionality in Eqs. ([M)) and (|55|). which is ln2 • 




FIG. 2: The estimated efficiency of the engine based on expected values of heats exchanged with the reservoirs (Eqs. (|38|) and 
(|39|l '). For near equifibrium, upto first order, the efficiency is 1/3 of Carnot value. 



10 

(ln(amax/amin)) , Can be related with heat capacity. The expected vahie of initial heat capacity of system i, defined 
as 

C, = / "^'^^ aU{ai)da„ (42) 

where we know 

Upon solving, the expected heat capacity in the initial state of the system is given exactly by 



C. 



, , ftp 

in 



m 



_T,(l + e'^»-/^0 T,(l + e°"-/^-) ^i + ga^ax/T. 

Then in the asymptotic limit, the leading term yields 



(44) 



— — In 2 

C, = C^ !"^ . . (45) 

This limiting value is independent of temperature of the system and thus indicates an analogy with a constant heat 
capacity thermodynamic system. 

Closing this section, we note that the requirement of consistency between the results of A and B implies, in 
an asymptotic limit, that the behavior expected from minimal prior information is the one which shows simple 
thermodynamic features such as constant heat capacity and equality of subsystem temperatures upon maximum work 
extraction. In the next section, we revisit the above analysis, but in much simplified form by considering an additional 
constraint. 

VI. ANALYSIS AT GIVEN THERMAL EFFICIENCY 

Let us impose the additional constraint by fixing the value of engine efficiency 77. Then there is essentially one energy 
scale say ai, to be specified in the model, because the other scale 02 is determined from the ratio 02/01 — (1 — ry). 
Let observer A treat oi as the uncertain parameter and denote the prior as 7r(ai). The second observer B chooses 02 
as the uncertain parameter and the corresponding prior as 7r*(a2). In other words, instead of ai, B parametrises the 
uncertain quantity as ai times some given constant which in the present case is (1 — r/). The probabilities assigned 
by A and B for a given choice of ai and 02, must satisfy 

7r(ai)dai = 7r*(a2)(ia2. (46) 

On the other hand, one should expect the same functional form for the prior in both cases, which is saying essentially 
that both A and B are in an identical state of knowledge, implying tt ~ tt* [3| . Thus relation ([IS)) is rewritten as 

7r(ai) = (1-77)^(01(1-77)), (47) 

a functional equation whose solution is given by ■n{x) oc 1/x. Thus with the additional constraint of a given efficiency, 
the prior can be assigned as the Jeffreys' prior. 
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Let us now illustrate the calculation for observer A. The appropriate normalised prior we have to consider is 
7r(ai) = [In (amax/amin)] (l/*^!)- From Eq. Q, the average work per cycle rewritten as function of ai and rj is given 
by 



W{ai,r]) = ail] 



1 



The expected work estimated by A is defined as VF(A) = J W{ai,ri)TT{ai)dai. After calculation we have, 



(48) 



W{A) = 



In 



1 -1 



Ta 



In 



1 + gQmax(l-'))/T2 



Ti In 



./Ti 



„/Ti 



(49) 



(1 - 77) '" V 1 + e'^""^^"''^/^" 

On the other hand, from the perspective of observer B who treats 02 as the unknown parameter, the corresponding 
prior is 7r(a2) = [In (amax/flmin)] (l/a2)- Upon writing the work per cycle as W{a2, rj), i.e. function of 02 and r/, we 
define W{B) ~ J W{a2,'i])TT{a2)da2, which is explicitly given by 



W{B) 



In 



V 



To 



■In 



1 



oa„„/T2 



Tiln 



I _^ ga„„x/(l-r;)Ti 



(50) 



In this case also, we see a difference in the average work expected by A and B, which is at variance with our assumption 
that A and B are in an equivalent state of knowledge. But again, we observe that in the aymptotic limit as considered 
in Section (fV)) . both the expressions for work get reduced to the following simpler form 

In 2 f^ T2 



W{A) K. W{B) 



V[Ti- 



(51) 



It is interesting to observe that the expected work in the asymptotic limit attains its optimal value at the well known 



Curzon-Ahlborn efficiency, r; = 1 — -y/Ta/Ti. This value, near equilibrium is approximated as 

(1-61) (1-6*^2 



V 



0(1- 



(52) 



2 8 

Note the difference in the above dependence on 9, with the one estimated by Eq. (|1T|) for the case of ignorance about 
two parameters. 

Finally, we may compare the amount of work expected in the general case (as calculated from Eqs. ((38)) and (jSQ]) ) 
when both oi and 02 are uncertain, with the special case when the efficiency is fixed a priori. In the general case, the 
expected efficiency of the engine is given by Eq. PO]) . so it is reasonable to take this value of efficiency in Eq. ([5T|) 
while making the above comparison. As shown in Fig. 3, the work in the general case is always less than the work in 
the special case and their ratio approaches the value of 3/4 as the temperature gradient goes to zero {9 -> 1). 



VII. CONCLUSIONS 



The main issue of this paper is about predictability in models of quantum thermodynamic machines in the presence 
of incomplete knowledge about the working medium. The implicit reason for uncertainty considered here is different 
from other scenarios usually considered: The source of this uncertainty does not lie in some unknown intrinsic dynamics 
of the system or due to some environment induced fiuctuations of the control parameters. We have not assumed any 
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FIG. 3: Expected work in the general case, is less than the expected work in the special case for which the efficiency is fixed 
a priori. Both expressions for work are calculated at the same efficiency as given by Eq. (I40|l . As S — s- 1, the ratio approaches 
a constant value 3/4. 



such objective mechanism which may presum.ably cause these values to become uncertain. Also, the uncertainty does 
not stem from lack of information about the preparation procedure for the quantum state of the working medium. 

We have assumed a Bayesian perspective about this (subjective) ignorance of the observer. In our model, there are 
two intrinsic energy scales of the working medium, which are under macroscopic control, for example, via external 
magnetic field. In our case, the hamiltonian itself is uncertain in as much as the values of control parameters entering 
in its description, are not pre-assigned. 

Once we assume a Bayesian perspective on our prior ignorance about the system, then the task is to narrow down the 
choice of prior for the uncertain parameters. The prior so assigned is supposed to quantify the given prior information 
such as the knowledge that these parameters are positive valued and the criterion for operartion as a heat engine. 
We make use of a consistency argument which says that different observers in possession of similar knowledge must 
assign similar priors. Invoking two observers A and B, an appropriate prior which we arrive at is the 1/x prior. Let 
us emphasize that in case further constraints or information from data become available, then the prior should be 
updated to include that additional information, by using suitable procedures 3j such as maximum entropy principle, 
or Bayes theorem. In the absence of data from observations, in our opinion, the initial prior so assigned has to be 
used to make inferences about physical quantities. The estimates for these quantities have been defined as the average 
value over the chosen prior. 

Further considerations lead us to investigate a particular asymptotic limit, because to maintain consistency, the 
observers A and B should arrive at similar estimates for a given quantity, if each is in an equivalent state of knowledge. 
It is in this limit, we observe classical thermodynamic features for the estimated quantities of our quantum heat engine. 
In particular, the expected mean energies of the two-level systems become proportional to their temperature, with the 
expected heat capacity C becoming independent of the temperature. It is also interesting to observe that the factor 
ln(aniax/amin) occuriug in the normalisation of the prior can be expressed in terms of the expected heat capacity of 
the system. The estimated efficiency behaves like one-third of Carnot value, for close to equilibrium (nearly equal 
bath temperatures). 

Further, if an additional constraint in the form of given value of efficiency is imposed, we observe that the number 
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of uncertain parameters in the problem reduces from two to one. It is then found that the expected work per cycle 
becomes optimal at CA efficiency. In Ref. 8], it was found that CA value is a lower bound for the efficiency at 
maximum work for this system. Proximity to this value at the optimal performance has been studied in many recent 
models and a certain universality has been observed [9|, Il0l |. In this context, this predictability about the optimal 
behavior by using an appropriate prior for ignorance, is remarkable. 

For the present model of heat engine, we have compared the estimated work per cycle in the general case of ignorance 
of two parameters with the special case of ignorance of one parameter only. Hcuristically, one can anticipate that the 
expected work should increase when additional constraints are introduced which lead to reduced uncertainty about 
the system. In particular, the ratio of these two values of work, tends to a limiting value of 3/4 as the two bath 
temperatures approach each other. 

Similarly, the estimated efficiency close to equilibrium for the general case is one-third of Carnot value, as contrasted 
with the one-half Carnot value as found for the special case. This may indicate that the general case of prior ignorance 
about two parameters may belong to a different universality class than the special case of ignorance about a single 



parameter, which is bounded from below by half Carnot value 10l.llll|. Further, we note that in a model of irreversible 



Brownian heat engine 



12l |. when the power is optimised with respect to the load and the barrier height, the efficiency 



at optimal power is found to be given by 

^ " 3-20(1 + In^) -02- i^^^i 

Close to equilibrium, this efficiency has the same expansion as given by Eq. (j4ip . upto second order. 

Concluding, our analysis suggests an unexpected relation between the average performance of certain infinite time 
quantum models of heat engines infered by a Bayesian analysis and the optimal performance of finite time thermo- 
dynamic models. In our opinion, the present study raises basic issues and open problems which we are addressing in 
future work, such as how robust is the predictability about efficiency at optimal work by choosing other models of 
engines, and to further investigate the influence of constraints on the expected behavior. 
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